In-plane dipole coupling anisotropy of a square ferromagnetic Heisenberg monolayer 
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In this study we calculate the dipole-coupling-induced quartic in-plane anisotropy of a square 
ferromagnetic Heisenberg monolayer. This anisotropy increases with an increasing temperature, 
reaching its maximum value close to the Curie temperature of the system. At T = the system is 
isotropic, besides a small remaining anisotropy due to the zero-point motion of quantum mechanical 
spins. The reason for the dipole-coupling-induced anisotropy is the disturbance of the square spin 
lattice due to thermal fluctuations ('order- by-disorder' effect). For usual ferromagnets its strength 
is small as compared to other anisotropic contributions, and decreases by application of an external 
magnetic field. The results are obtained from a Heisenberg Hamiltonian by application of a mean 
field approach for a spin cluster, as well as from a many-body Green's function theory within the 
Tyablikov-decoupling (RPA). 

PACS numbers: 75.70.-i, 75.50.Ee, 75.70.Cn 
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I. INTRODUCTION 

The investigation of the magnetic properties of fer- 
romagnetic ultrathin films is a field of intense current 
interest!! In this contribution we determine the mag- 
netic in-plane anisotropy of a square two-dimensional 
(2D) Heisenberg ferromagnet in presence of the magnetic 
dipole interaction. In 2D magnets the action of magnetic 
anisotropics is-,twpfold. First, they induce a long range 
magnetizationcltJa with a Curie temperature of the order 
of the exchange coupling^, whereas a strictly isotropic 
2D Heisenberg magnet does not exhibit any long range 
magnetiCj-prder at finite temperatures (Mermin-Wagner- 
theoremp. Moreover, the long range character of the 
dipole interaction itself is sufficient to stabilize the mag- 
netization in a 2D magnet. Secondly, the anisotropics 
determine the easy and hard axes of the magnetization 
with respect to the lattice framed. For a ferromagnetic 
thin film the dipole coupling prefers an in-plane direction 
of the magnetization (shape anisotropy, demagnetizing 
field). Except for very small thicknesses thin films are 
usually in-plane magnetized. The corresponding in-plane 
easy axis depends on the lattice symmetry of the film 
face. For a rectangular (110) monolayer the dipole in- 
teraction induces a uniaxial in-plane anisotropy, the easy 
axis being the in-plane axial direction with the smaller 
lattice constant, as can easily be shown by calculating 
the corresponding lattice sums. 

On the other hand, for a square monolayer the dipole 
coupling, although itself not rotational invariant, ex- 
hibits a continuous energy degeneracy for classical spins 
at T = 0. Thus, due to the Mermin- Wagner theorem a 
long range magnetic order is not expected at finite tem- 
peratures. However, as has been shown by Monte Carlo 
calculations and interacting spin wave theory, a magnetic 
ordering and a critical temperature exists for dipole cou- 
pled spins on a square latticeQ, since the magnetic exci- 
tations are not continuously degenerate. In this case a 
quartic in-plane anisotropy is present, the corresponding 
easy axes being the edges of the square lattice. In other 



nomenon is an exam 
in frustrated magnetf 



words, the density of states and thus the entropy depends 
on the magnetic direction within the lattice. This phe- 
e of the 'order-by-disorder' effect 
. In these systems thermal fluctu- 
ations or lattice disturbances partly remove frustrations, 
and a collective magnetic ordering emerges at finite tem- 
peratures. 

Whereas the existence of a finite magnetization and 
a quartic anisotropy at finite temperatures for dipole 
coupled spins on a square lattice has been proven, the 
strength and the temperature behavior of this anisotropy 
has not been determined yet. Of particular interest is 
whether the dipole-coupling-induced quartic anisotropy 
can be measured, and how this anisotropy compares with 
the single-ion quartic lattice anisotropy due to the spin- 
orbit coupling. In this study we will adress these ques- 
tions with the help of a mean field approximation, taking 
into account a square spin cluster with different num- 
bers of spins (Oguchi-theory)Bt2l, as outlined in Sec. II. 
Furthermore, we apply a many-body Green's function 
theory within the Tyablikov-decoupling (RPA) by con- 
sidering collective magnetic-eptcitations (spin waves) and 
interactions between themUO, which are in particular 
important for 2D magnetic systems. Our results are sum- 
marized and discussed in Sec. III. 



II. THEORY 



The free energy F(T, (f>) is determined within a 
Heisenberg-type Hamiltonian as a function of the tem- 
perature T and the in-plane angle <j). A square (001) 
ferromagnetic monolayer is considered, which is spanned 
by the xz- plane, the z- direction refers to <fi = 0. The 
isotropic nearest neighbor Heisenberg exchange coupling 
J, the magnetic dipole-dipole interaction, and an exter- 
nal magnetic field B between localized Heisenberg spins 
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Si on lattice sites i are taken into account: 



j - g^B 



E BS * 



+^r ! E k K s * s ^ - 3 ( r « Si ) ( r « s ^)l ' (1) 



with g the Lande factor and fis the Bohr magneton. The 
spin quantum number is assumed to be 5 = 1/2. The 
lattice vector between sites i and j is given by ry, with 
r ij = l r ij'l the distance. All spins are assumed to be 
aligned parallely (ferromagnetic order), the magnetiza- 
tion m;(T) = m(T) = (0,0, m(T)) is directed along the 
magnetic field B = (0, 0, B). 

To account for a varying angle of the magnetization 
with respect to the lattice frame, the lattice is rotated by 
the in-plane angle <p with respect to the z- axis, yielding 
the rotated lattice vectors: 



r-j sin< 



r-j cos< 



rfi cos < 



(2) 



Clearly, the only source in Eq.(|l|) for a possible in-plane 
anisotropy is the double scalar product of the dipole 
interaction. Within thej-single-spin (Bragg- Williams-) 
mean field approximationtHI the dipole coupling yields an 
additional contribution to the molecular field. The re- 
sulting single-spin Hamiltonian reads 



H (1) = — (qJm(T) + w 5(0, 0) m(T) + g/j, B b\ S* 

+jm 2 (T)(gJ + W S(0,0)), (3) 

with q — 4 the number of nearest neighbors (coordina- 
tion number), w = (y/is) 2 /ciq the strength of the dipole 
interaction, a the lattice constant, 5(0, 0) ~ 4.51681 the 
corresponding lattice surrJlj, and m(T) = (Sf ) the mag- 
netization. Thus, within this approximation the dipole 
interaction is isotropic for a square monolayer and does 
not depend on the in-plane angle <f>. 

Our approach will nowJae improved by applying the 
so-called Oguchi methocH'El A number N > 1 of neigh- 
boring spins is considered, the interactions between the 
N spins in this cluster are treated exactly. The remain- 
ing lattice is coupled to the spin cluster by a molecular 
field. To conserve the symmetry of the square lattice, 
we take into account only square-shaped spin clusters, 
the three smallest possible clusters are characterized by 
N = 4, 9, and 16. Note that the resulting cluster Hamil- 
tonians HS N > are usually non-diagonal, i.e. contributions 
proportional to the S x - and S y - spin components may be 
present. Therefore, as a further approximation we con- 
sider only the diagonal elements. This approach guaran- 
tees that the magnetization is directed along the external 
magnetic field. As an example, the Hamiltonian for the 
four-spin cluster (N — 4) reads: 



H (4) = -( J - w) (5 1 2 5 2 2 + 5f SI + S z 2 Sl + 5 3 2 5f ) - ^(Sf 5| + 5f 5f) 

-3 w sin 2 (S{S Z 2 + 5|5|) + cos 2 cj, (SfSf + S z 2 Sl) + 2 {l 

- (2 J m(T) + g[iB B) + 5 2 z + 5f + 5|) - w m{T) (S z + 5|) (5(0, 0) 
-wm(T) {S z 2 + SI) (5(0,0) - 1 - -6 cos <t> sin^)) +4 Jm 2 (T) + 2wm 2 (T) (5(0,0) - 1 - . (4) 



1 ~ 25/2^ + 6 cos< ^ sin< 



r 



The partition function for this case is given by: 

+s 

Z^{T,B,4>)= E exp(-/?HW), (5) 

S 1 ^ , S| , S3 , S| — — s 

and the (average) magnetization by: 

1 



m(T, B, 



AZW{T,B, 



+s 



E (S? + Si + SI + SI) exp(-(3n {4) ) , (6) 

Sf ,S% ,Sf,3^=— S 

with [3 — 1/ksT and fee the Boltzmann constant. The 



corresponding expressions for the other investigated sizes 
are rather lengthy and not shown here. The difference 
of the free energy F(T,B,<f>) = -k B T lnZ(T,B,(j>) be- 
tween the diagonal (</> = tt/4) and the axial directions 
((f) = 0) yields the in-plane anisotropy K,^dip(T, B), which 
will be calculated as a function of the temperature T and 
the external magnetic field B. 

It is important to mention that the mean field theory 
as described in the preceeding subsection does not fulfill 
the Mermin- Wagner-theorem for the isotropic 2D Heisen- 
berg magnet. In this case the long range magnetic order 
becomes unstable againstxollective magnetic excitations 
with long wavelengthdSffl. Thus, the consideration of 



these spin waves is very important for the magnetic be- 
havior of a ferromagnetic monolayer. For the determina- 
tion of the dipole-coupling-induced in-plane anisotropy 
we apply in addition a many-body Green's function the- 
ory to the Hamiltonian Eq. (fil) , after the lattice is rotated 
by the in-plane angle <fi, Eq.l||) . We consider the Green's 



functions Gf A 



which are 



solved within the Tyablikov-decouplingjlil. Since interac- 
tions between spin waves are partly taken into account, 
the magnetic properties can be determined up to the 
Curie temperature by this method. The respective for- 
malism is described in detail inB, thus we outline merely 
necessary extensions in the Appendix. The magnetiza- 
tion m(T, B, <j>), the internal energy E(T, B, (j>), and the 
free energy F(T, B, (j>) are calculated. For comparison, 
we determine the corresponding quantities alse-by con- 
sidering the Holstein-Primakoff approximationE-3, which 
is valid at low temperatures. 



III. RESULTS AND DISCUSSION 

By application of the described methods we calcu- 
late the effective dipole-coupling-induced in-plane mag- 
netic anisotropy IC^.dipiT, B) for a square ferromag- 
netic monolayer as a function of the temperature T 
and the applied magnetic field B. A finite value in- 
dicates that this anisotropy is caused by the mag- 
netic dipole coupling, since other sources of magnetic 
anisotropies are not present in the Hamiltonian Eq. ([!]). 
For the strength of the dipole interaction we choose 
(5Mb) 2 / o = w = 0.01 J, which is appropriate for 3d- 
transition metal ferromagnets. JC^dip(T, B) is given in 
units of the energy difference between the perpendicular 
and the in-plane magnetization (demagnetizing energy) 
E dema g(0) = (3/2)uiS'(0,0)m 2 (0) k2ttwS 2 for a ferro- 
magnetic monolayer at T = 0. 

The free energy F(T, <fi) exhibits a four-fold symmetry 
as a function of the in-plane angle <$>. The easy mag- 
netic axes are directed along the edges (4> = 0, it/2, etc.), 
and the hard axes along the diagonals = tt/4, 37t/4, 
etc.) of the square lattice as obtained inQ. The quartic 
anisotropy JC^dipiT, B) is depicted in Fig.l as a func- 
tion of the temperature, calculated by the Oguchi ap- 
proach. Three different sizes of the spin cluster are con- 
sidered (N = 4, 9, 16), as well as different magnetic field 
strengths B. The Curie temperature Tc{N) decreases 
with an increasing N, since additional spin correlations 
are considered. As can be seen, K,^. d ip(T, B) increases 
with an increasing temperature, reaching a maximum at 
Tc- This temperature behavior is in striking contrast to 
the corresponding behavior of other anisotropic contri- 
butions. Usually the effective anisotropies decrease with 
an increasing temperature, and vanish at Tc for B = 0. 
l^4,dip(T, B) exhibits a cusp at Tc and decreases for tem- 
peratures T > Tc- By increasing the size N of the spin 
cluster, K±.dip{T, B) becomes larger. An application of 
a magnetic field reduces IC^dipiT, B), and the cusp near 
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FIG. 1: Effective dipole-coupling-induced quartic in-plane 
magnetic anisotropy /C4, dip (T, B) for a square monolayer 
as a function of the temperature T, calculated within the 
Oguchi mean field method for spin cluster sizes TV = 
4, 9, 16. Two different magnetic field strengths B are con- 
sidered (gusB /E dEmag (0) — 0, 1.5). The field energy and 
K,4 t dip(T, B) are given in units of the demagnetizing energy 
Edemag(0)- The dot-dashed line shows Ki,dip(T, B) for B = 
and N = 4 by considering non-diagonal matrix elements of 
the Hamiltonian, cf. Sec. II. For comparison, also the effective 
single-ion quartic anisotropy K.4 iSO (T) resulting from the spin- 
orbit interaction is shown, calculated within a thermodynamic 
perturbation theory by assuming £4,30(0) = 0.01 E de mag(0) 
(dotted line). 



Tc changes to a maximum. For comparison, in Fig.l 
we show /C4,di P (T ', B) for B = and N — 4 by consid- 
ering also non-diagonal elements in the cluster Hamilto- 
nian H.( 4 \ cf. Sec. II. A decrease of Tc and an increase 
of K,4 : di P (T, B) by ~ 10% is obtained with respect to the 
calculations which consider diagonal elements only. 

In addition, we have calculated IC^dip{T, B) by ap- 
plication of_ the many-body Green's function theory as 
outlined ino and in the Appendix. For w = 0.01 J 
the Curie temperature is calculated to be kBT^ PA — 
0.373 J. In Fig.2(a,b) we present the resulting magneti- 
zation m(T,B) and anisotropy K,4 td i P (T, B) as functions 
of the temperature and the magnetic field. In accor- 
dance with the results obtained from the Oguchi ap- 
proach, K-4^dip(T, B) increases with an increasing tem- 
perature. Due to the consideration of quantum me- 
chanical spins a finite value of K,^dip{T,B) is already 
present at T — 0, resulting from the zero-point spin mo- 
tion (quantum fluctuations). Note that the presence of 
the dipole coupling causes a-,pen-saturated magnetiza- 
tion m(T = 0) < S = l/2tM In Fig.2 also the re- 
sults as obtained from the Holstein-Primakoff approxima- 
tion are shown, which neglects interactions between spin 
waves. The corresponding magnetization m(T, B) and 
anisotropy K.4 td i P (T, B) at low temperatures are close to 
the results as obtained from the RPA. At elevated tem- 
peratures T > Tc/3 the Holstein-Primakoff approxima- 
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FIG. 2: Magnetization m(T, B) (a) and effective quartic 
in-plane magnetic anisotropy K,4 t di P {T, B) (b) for a square 
monolayer as functions of the temperature T and the mag- 
netic field B, calculated by application of the Green's function 
method within the Tyablikov-decoupling (RPA). The Curie 
temperature is calculated to be ksTQ PA = 0.373 J. The 
field energy and K,i : dip{T, B) are given in units of the de- 
magnetizing energy Edemag{0)- In addition, the results as 
obtained from the Holstein-Primakoff approximation (HP) for 
5 = are depicted (dashed lines). For comparison, the dotted 
line shows the effective single-ion quartic anisotropy K,4, t30 {T), 
with /C 4 , so (0) = 0.01 E demag (0). 

tion is no longer valid. 

As has already been obtained from Moiaie Carlo cal- 
culations and interacting spin wave theoryu, the dipole- 
coupling-induced in-plane anisotropy is caused by mag- 
netic fluctuations. This order-by-disorder effect is. thus 
controlled by the magnetic entropy of the systemtl. An 
energetic influence should become apparent at T — al- 
ready within a mean field approximation. Besides the 
effects resulting from the zero-point motion, the terms 
of the dipole interaction dependent on the in-plane an- 
gle cf> cancel for a strictly square magnetic lattice, i.e. 
for a fully saturated ferromagnetic state. At finite tem- 
peratures thermal fluctuations disturb this square peri- 
odicity, causing an effective quartic in-plane magnetic 
anisotropy, which reflects the underlying spatial lattice 
symmetry. The single-spin mean field theory conserves 
the periodicity of the square magnetic lattice, since all 
other spin operators are replaced by uniform expecta- 
tion values. Thus, within this approximation one obtains 
JCi t dip(T, B) = for all temperatures. On the other hand, 
within the Oguchi approach the interactions between a 
number of spins are treated exactly, considering a few col- 
lective magnetic excitations with short wavelengths rang- 
ing over several lattice constants, and resulting in a finite 
value for K^ap^T, B) for T > 0. The many-body Green's 



function theory takes into account spin waves with all 
possible wavelengths. Since K,^dip{T,B) is maximal at 
elevated temperatures, K^dip{T, B) is caused by mag- 
netic fluctuations in particular with short wavelengths, 
which become excited in this temperature range. 

Within the Oguchi approach we calculate a maximum 
value of JC^dipiT, B) of the order of 0.1% of the demagne- 
tizing energy E dem ag{0)- On the other hand, }C^ d i P (T, B) 
as calculated from the Green's function method is 
more than a magnitude larger, namely K,^dip(T, B) ~ 
0.01 Edemag(0) near Tq- Similarly, the action of a mag- 
netic field B on the anisotropy IC^dipiT, B) is calculated 
to be much more pronounced within the Green's func- 
tion theory than within the Oguchi approach, cf. Figs.l 
and 2. The reason for the strong differences between the 
results obtained from the two methods is that the former 
one considers collective magnetic excitations with long 
wavelengths. These are known to have slmng effects on 
the magnetic properties of ultrathin filmscm. In this case 
the magnetic field acts merely on spin blocks, i.e. cor- 
related regions of neighboring spins characterized by the 
short range order parameter!!!], resulting in a much larger 
magnetic response. For instance, the induced magnetiza- 
tion in a ferromagnetic trilayer has been determined to 
be anpOrder of magnitude larger by consideration of spin 
waveslij. The collective magnetic excitations are most 
pronounced for a single magnetic layertl, the zero-point 
motion is strongest for the spin number S = 1/2. Note 
that the theoretical methods applied in the present study 
are much less demanding than Monte Carlo simulationsQ. 

The obtained dipole-coupling-induced in-plane 
anisotropy is small as compared to other (effective) 
anisotropics, since we have considered interaction 
strengths appropriate for 3d- transition metal ferromag- 
ncis. A small value for K-4 t dip(T, B) has been conjectured 
irJ3. If a ten times larger ratio w/J between the dipole 
coupling strength and the exchange interaction is 
assumed, the ratio K,A,dip{T, B)/ Ed e mag(0) increases 
roughly by the same factor. Since usually the effective 
anisotropics decrease with an increasing temperature, 
K-A,dip{T,B) might become observable at elevated tem- 
peratures. For comparison, we show in Figs.l and 2 
also the effective single-ion quartic in-plane anisotropy 
fci.so(T) resulting from the spin-orbit interaction. This 
quantity is calculated, with the help of a thermodynamic 
perturbation theoryO. By assuming its strength to be 
^Q.so(O) = 0.01 Edemag(O), and adapting corresponding 
Curie temperatures, lCi tSO {T) becomes comparable to 
fC4,di P (T,B) at T/T c ~ 0.8 - 0.9. In this temperature 
range the total quartic in-plane anisotropy should 
increase again, or exhibit an in-plane magnetic reori- 
entation. As mentioned, this dipole-coupling-induced 
anisotropy should be more apparent for ferromagnetic 
(001) thin films with a large w/J- ratio. 

We like to comment on the fact that a finite anisotropy 
is obtained for T > Tq. Usually the effective anisotropy 
as observed for a collectively ordered ferromagnetic 
state vanishes above Tq- However, a vanishing effec- 
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tive anisotropy for T > Tc does not indicate that the 
anisotropy as present in the Heisenberg Hamiltonian dis- 
appears. A single (paramagnetic) spin is still subject to 
e.g. the single-ion uniaxial anisotropy Ki also if the net 
magnetization is zero. The resulting free energy differ- 
ence between the easy and the hard magnetic directions 
('paramagnetic anisotropy') behaves as oc K^/ksT for 
K2 <C ksT, and is rather small if K2 is small as compared 
to the exchange interaction J oc ksTc- With regard to 
the present study, a finite value of K.4 t ,a P (T,B — 0) for 
T > T c reflects, the dipole-coupling- induced anisotropy 
of a spin blockta. 

Note that the free energy F(T, B) as obtained from 
the Green's function method yields unphysical results for 
temperatures near and above Tc- As can be shownlia, 
for large temperatures the free energy F(T, B) as cal- 
culated by this method does not approach the value 
— fcsTln2, which is the free energy of a single, non- 
interacting S — 1/2- spin. This unphysical behavior re- 
sults in a still increasing K.^dip{T,B) for T > Tc, cf. 
Fig. 2(b). Nevertheless, we expect that the Green's func- 
tion approach reflects the correct behavior of K.4, t dip(T, B) 
for T < Tc, since its behavior is corroborated by the 
Oguchi method and the Holstein-Primakoff approach. 

We have considered a square monolayer only. A similar 
effect is expected also for thicker films with a square (001) 
face, or for a hexagonal (111) thin film. In addition, a 
corresponding effect may emerge for a three-dimensional 
cubic lattice, for which the dipole interaction cancels 
exactly. At elevated temperatures, however, caused by 
thermal fluctuations the frustration due to this period- 
icity will be lifted. The dipole coupling then induces a 
cubic magnetic anisotropy with easy axes directed along 
e.g. the edges of the cubic lattice. Also is this case an in- 
creasing anisotropy with an increasing temperature is ex- 
pected. To obtain an anisotropic contribution caused by 
the magnetic dipole interaction via the order-by-disorder 
effect, we emphasize that the lattice symmetry must be 
strictly cubic or square. A lattice distortion e.g. due to 
magnetostrictive effects will probably blur this effect. 

P.J.J, gratefully acknowledges the invitation and the 
hospitality of the I. Institute for Theoretical Physics of 
the Hamburg University. Fruitful discussions with Prof. 
K.D. Schotte, Free University of Berlin, are acknowl- 
edged. 



Appendix: Many-body Green's function theory 



The following commutator Green's functions in the 
frequency space are taken into account: Gfj(u>) — 
{(Sf ; S~)) u , a = +,—,z. By considering the Hamil- 
tonian Eq. ([[]), these Green's functions are— salved in 
the usual way by the equation of motionaOO The 
higher order Green's functionS|— are approximated by 
the Tyablikov-decoupling (EPA)!!!!, e.g. ({S%Sf;Sr}) ~ 
(S Z )((S^~; S~)), resulting in a system of linear equations 



for the Gfj(uj). By application of a Fourier transforma- 
tion into the 2D wave vector space k|| = k, one obtains 
G z {k,uj) = 0, if (S+) = (S~) = 0. The remaining two 
Green's functions G ± (k, uj) are easily evaluatec 



G + (k,w) = 2m(T) 
G-(k,uj) = 2m(T) — 



to + a(k) 



w 2 -e 2 (k) 
2 (k) 



(7) 
(8) 



with m(T) = (S z ). The magnon dispersion relation e(k) 
is given by 



e 2 (k) = a 2 (k)-6 2 (k), 
a(k) = gfi B B + m(T)(j(0)- J(k) 



(9) 
(10) 



- w m(T) 



S(0,0) + S(k z ,k x ){l-^cos 2 0) 



+ S(k x ,k z ) (l - ^sin 2 0) +3T(k z ,k x ) cos0 sin</> 



b(k) = ^wto(T) S{k z ,k x ) cos 2 + S{k x ,k z ) sin 2 



2 T(k z , k x ) cos 4> sin i 



(11) 



J(k) — 2 J (cos k x ao + cos k z ao) is the Fourier transform 
of the exchange interaction,r-with a the lattice constant. 



13; 



The oscillating lattice sumstia are defined by 



S(k x , k z ) 
T(k z , k x ) 



E 

1— — 
+00 

E (12 + n 2)5/2 



l,n— — 00 
+00 



[l 2 + 71 2 ) 5 /2 



,n— — 00 



exp(—ik x aol — ik z aon) , 

exp(—ik z aol — ik x aon) , 
(12) 



where the terms with I = n = have to be omitted. 
The following expressions are valid for the spin quantum 
number S = 1/2. The magnetization m(T) is obtained 
from the spectral theoremO: 

m(TA) = I-l^(5-5+)(k) 



N 

1 m(T) 

2 W 



E(fi oath* -l), (13) 



v 6 (k) 



with x = e(k)/2fcsT. The summation extends over the 
first Brillouin zone, N denotes the number of k-points. 
By calculating the expectation value (S Z S Z ) with the 
help of G + (k, u> we obtain for the internal energy 
E(T, 4>) per spin: 



E(T,4) = 

s(k) 



"(e 2 (k) 



f e(k)) cothrr - a(k) + 6(k) 



coth x — 1 



(14) 
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with the denotations 

Eo = ~gn B B-^(j{0)+wS(0,0)) , (15) 
ei(k) = glMB B + ±(j(0)- J(k)) 

+| (5(0, 0) + ST(fc„ fe x ) + S(fc x , fc,)) , (16) 
e 2 (k) = ij(k)-^(s(k z ,k x ) + S(k x ,k z )y (17) 



Finally, the free energy F(T, <fkis calculated from a tem- 
perature integral over E(T, </>)E3. The Holstein-Primakoff 
approximation!!-! is obtained by replacing the magnetiza- 
tion m(T) on the right side of Eq. (|l3|) by its saturation 
value 5 = 1/2. 
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